Fermi ID quantum gas: Luttinger liquid approach and spin-charge separation 
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We discuss the properties of quasi- ID quantum gases of 
fermionic atoms using the Luttinger liquid theory, includ- 
ing the presence of an optical lattice and of a longitudinal 
trapping potential. We analyze in particular the nature and 
manifestations of spin-charge separation, where in the case 
of atoms "spin" and "charge" refers to two internal atomic 
states and the atomic mass density, respectively. 



I. INTRODUCTION 

The recent experimental possibilities to trap ultra-cold 
quantum gases in quasi-lD magnetic and optical wave- 
guide aroused a renewed interest in quantum liquids with 
reduced dimensionality and their peculiar properties. In 
this context, most of the recent theoretical and exper- 
imental work has focused on quasi- ID Bose gases and 
Bose-Einstein condensates [1,2]. In particular, the equi- 
librium properties have been already discussed in several 
papers (see Ref. [3] and references therein), the phase 
fluctuations of the classical field (order parameter) de- 
scribing a quasi-lD condensate have been predicted in 
Ref. [4] , and experimental evidence was reported in Ref. 
[5]. In Ref. [6] the collective excitations for ID Bose gases 
have been calculated for different configurations, ranging 
from the ID mean-field to the Tonks regime. 

At the same time also the possibility to include an 
optical lattice has become available, opening new per- 
spectives in the study of coherence phenomena [7] . For 
instance Bloch oscillation have been seen in Ref. [8]. In 
Ref. [9], it is claimed that a direct measurement of the 
critical Josephson current has been obtained, and very re- 
cently the intriguing quantum phase transition between 
superfluid and Mott-insulator, predicted some years ago 
in Ref. [10], was observed for a Bose gas loaded in a 
3D optical lattice in Ref. [11]. Also the dynamics of a 
trapped Bose-Einstein condensate in the presence of 2D 
and ID optical lattice has been calculated in Ref. [12]. 

On the other hand remarkable progress in trapping 
and cooling Fermi gases, has been made to reach the de- 
generate regime [13]. On the theoretical side different 
aspects of (harmonically) trapped Fermi gases have been 
explored, often in analogy with what has been done for 
Bose gases [14-16]. Experimentally, measurements of the 
collective modes of degenerate Fermi gases have been re- 
ported (see Ref. [17] and reference therein). 

Fermi cold-gases, besides having an interest per se, can 
represent a useful tool to implement solid state systems 
and to test fundamental theories, often grown and de- 



veloped in a different area. Trapped gases are cleaner 
and more flexible systems than the solid state counter- 
part, where the experimentally most relevant ID sys- 
tems are semiconductor quantum wires and nanotubes. 
In addition, it is also possible to use magnetic flelds 
to largely modify the scattering properties of atoms in 
atomic clouds [18]. 

Returning to ID systems, what is interesting is that 
due to the reduced dimensionality the interaction, al- 
though small, strongly modifles the properties of the 
system. Low dimensions amplify the role of quantum 
fluctuations and enhance correlations. Both the ground 
state and the excitations exhibit strong correlation effects 
and posses a number of exotic properties, e.g., fractional 
statistics (see Ref. [19-21] and references therein). An- 
other fundamental interest in ID systems originates from 
the fact that there is a fair number of exactly solvable 
models, allowing the research of non-perturbative effects 
[22]. In particular Fermi-Landau theory which describes 
3D interacting fermions is not longer applicable (e.g., see 
[20]). Instead, in a quite wide range of parameters, ID 
quantum liquids can be described using the Tomonaga- 
Luttinger model [23], and in this case they are named 
Luttinger liquids, in contrast to (normal) Fermi liquids. 
The term Luttinger liquid was coined by Haldane [19] 
to describe the universal low-energy properties of gapless 
ID quantum systems. The Luttinger liquid language has 
been already applied to study the properties of (quasi)- 
ID Bose gases [24,25]. The goal of this paper is to study 
trapped ultra-cold Fermi gases in quasi-lD geometry us- 
ing Luttinger theory and to propose an experimental way 
to see a typical Luttinger liquid behaviour, the so-called 
spin-charge separation. Indeed one of the hallmarks of 
a spin-1/2 Luttinger liquid involves the complete separa- 
tion in the dynamics of spin and charge degrees of free- 
dom. Both branches, which exhaust (in a infinite system) 
the spectrum, are sound-like with different propagation 
velocities. This phenomena has never seen in a clear way 
in actual condensed matter systems (see e.g. [26]). 

The paper, which can be considered a self-consistent, 
extended version of Ref. [27] , is structured as follows. In 
SEC. II we introduce the models. Using standard argu- 
ment, we write down a bosonized Hamiltonian to describe 
the low-energy excitations of the systems (SEC. III). We 
used the available exact solutions to relate the so-called 
Luttinger parameters to the microscopic parameters of 
the system for small and large interaction (SEC. III). In 
SEC. IV we shall propose different experimental ways to 
see clearly spin-charge separation in our systems, and in 
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SEC. V we study the hfe time of the excitations. In SEC. 
VI we propose an implemetation of spinless Luttinger hq- 
uid using strongly interacting bosons in an optical lattice. 



II. PHYSICAL MODEL 

We consider a dilute gas of fermionic atoms of mass 
TO with two ground states, |(t =|, I), representing a spin- 
1/2, e.g. the states |F = 1/2, M^^ = ±1/2) of ^Li [28]. In 
the following we assume the number of atom in each level 
is the same: N^r = N/2. The atoms are confined in a very 
elongated trap along the a;-direction and cooled below 
the Fermi temperature Tp- The transverse confinement 
is considered to be harmonic with frequency lu±, strong 
enough so that the transverse degrees of freedom frozen. 
In order to have a (quasi) ID system the tight transverse 
trapping must exceed the characteristic energy scale of 
the longitudinal motion. Due to the quantum degeneracy 
the longitudinal motion has all the energy levels up to the 
Fermi-energy e^? ^ ksTp occupied, thus the requirement 
is ep 'ti huj±. 

At low temperature, due to the Pauli exclusion prin- 
ciple, only inter-species s-wave collisions are allowed. 
All the relevant interactions are thus characterized by a 
single parameter, the inter-component scattering length 
a. Introducing the harmonic oscillator transverse length 
l± — y^h/mLL)±, the effective ID interaction can be rep- 
resented as a zero-range potential with strength g — 
2TTli^a/ml\, when a ^ /j^ [2]. Thus the system can be 
described by the following ID Hamiltonian 
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(1) 



where ipa is the ID field operator for atoms in the state 
a. The external potential V^xt.o-, includes both the lon- 
gitudinal confinement Vl,^ and a possible superimposed 
optical lattice Vopt,a-- The longitudinal confinement is ei- 
ther a box of length L (we will also refer to this as the 
homogeneous case), or an harmonic potential (inhomo- 
geneous gas) with frequency co^, i.e. V^.o- = l/2TOw^a;^. 
The confinement potential can act differently on the two 
species, and this feature can be used in an experiment, as 
we will show below, to test the "spin-charge separation" . 
The simplest way to create a ID periodic potential for 
neutral atoms is by superimposing two linearly polarized 
counter-propagating laser beams. The created conserva- 
tive potential is Vopt,a = Vq sin^(fca;), FIG. (1), where the 
lattice period d — n/k — A/2 is fixed by the wavelength, 
A, of the laser light (see, e.g., [29]). We can expand the 
field operators in the first band Wannier basis 



w(x — Xi)Ci, 



(2) 



where wix — Xi) is the Wannier function centered at the 
i-ih site and Ci^a the annihilation operators for a fermions 
in the i-th site with spin a. Plugging the previous expres- 
sion into Eq. (1) and retaining at most nearest neighboor 
terms for the kinetic energy part we reduce the starting 
Hamiltonian to the Fermi-Hubbard Hamiltonian 

HpH = - J ^(c|+i,aCi,cr + h.c.) + ^ ei,c,ni,c 

i.a i,(r 

+^X!"'»'T"-»,i' (3) 



counts the number of a- 



where operator Tii^a — a^t,a 
atoms in the site i. The parameter U correspond to the 
strength of the on site inter-species interaction depend- 
ing on J is the hopping term between adjacent sites 
and ei^a — VL,a{xi) describes an energy offset of each lat- 
tice site. Note that for d — > one recovers a continuum 
model. 
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FIG. 1. Schematic representation of fermions in an optical 
lattice. The parameter J is the amplitude to jump from one 
site to the adjacent one. U is the on-site energy due to the 
interaction between atoms in different internal levels. 



III. LUTTINGER LIQUID APPROACH 

In what follows we will often use the words charge and 
spin though in our system we do not have, strictly speak- 
ing, charge nor spin. We used them to maintain the jar- 
gon of the Luttinger physics. We identify the charge with 
the mass and the spin with the two atomic internal levels. 
Thus when we talk about charge density fiuctuations we 
mean density fluctuations, when we talk about spin den- 
sity fluctuations we mean relative density fluctuations, 
the so called composite excitations. 



A. Homogeneous case 

In this subsection we will review briefly what it means 
to describe the low energy properties of the system using 
the so-called Tomonaga-Luttinger model. 
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We consider the honiogcncous case, i.e., Vl.ct = 0, and 
we will study later the inhomogeneous case. The stan- 
dard way to represent the Eq. (1) for a gas in a box of 
length L with periodic boundary condition, without the 
underlying optical lattice, in the momentum representa- 
tion is 



H : 



'k,(T 



Ck,, 



E - 

ki,k2,q 



(4) 



where Sk = h^k'^/2m — /x. 

By means of renormalization group (RG), it is possi- 
ble to show that the low-energy properties of the sys- 
tem can be described using an exact solvable model, the 
Tomonaga-Luttinger model [19-21,23]. First of all in a 
ID system the Fermi "surface" is composed by just two 
points, zbfci?. In the non-interacting case one has ikp = 
±Trhn/2, where n = {N^ + Ni)/L is the total density of 
the gas. Let us, then, introduce the the left- (around 
—kp) and right- (around +kp) moving fermion annihi- 
lation (creation) operators Cr^aik), where r = R, L, and 
the respective densities Pr,s{<l) = X^fecJ.^(fc -I- q)cr,a{k)- 
In the model the four particle species {R 1 , R .[. L 1^ . L l) 
have unbounded free dispersion relation er,q = arhvpQ, 
where vp = hTm/2m is the Fermi velocity and clik^l) = 
-|-(— ); all the states up to the Fermi energy, in the same 
spirit of the Dirac sea, are occupied. Thus one has two 
branches of particles with boundless energy and momen- 
tum and the densities obey Bose type commutation rela- 
tion [20]. The model Hamiltonian written in the notation 
of the so-called "g-ology" (see for instance [30]), in real 
space representation, takes the form 

HlM = E / dXTrhVp{pR^a{xf +PL,a{^f) 

+ 54.1 E / dxpr,a{x)pr-a{x) 

/ dxi:\^^i)R,afl}^n_^i)L-a- (5) 



+ ; 



Specific names arc given to the parameters which iden- 
tify the strength of a particular scattering process. For 
instance g2_L is the strength of the forward scattering 
between particles belonging to different branches and 
with different spin state, {Ra,L — a) {Ra,L — a), 
while gi± is the strength of the backward scattering 
{Ra,L — a) ^ {La,R — a). In a complete model also 
"longitudinal" scattering {r(T,r'a) {ra,r'a) could be 
present. In our case, where only inter-species collisions 
are allowed, we have 52 1| =5411 = Cl- 
in the sense of RG in the long- wavelength limit, it is 
often the case that gi± is renormalized to zero [20,31-33]. 



This is true in our case. The final step to solve the prob- 



lem is to introduce four new boson fields 



n,,, with 



= c, s, such that 

Pr,cr{x) = 



[dx<Pc 



rUc + (T{d:r<ps - arils)]. (6) 



The bosonic fields are related to the fluctuations of the 
total density Pc = dx(l>c/V^, of the spin-density Ps = 
dx'Ps/V^i of the current density jc = —Tlc/^/^T and of 
the spin current density js = —Us/y/n. The fields <j) and 
n are conjugated and they satisfies the bosonic commuta- 
tion relation [(j>,y{x), (pi,/ (x')] ~ [n,y(.x), 11^' (.-r')] = and 
[(j)^{x),Il^i{x')] = iS^yS{x — x'). Substituting Eq. (6) in 
Eq. (5) one obtains the bosonized Hamiltonian [20] 



H 



dx 



(7) 



We are left with the Hamiltonian of two independent elas- 
tic strings with the eigcnmodes corresponding to the col- 
lective density and spin-density fluctuation of the fermion 
liquid. The parameters u^'s are the "sound" velocities, 
while the K^s are related to the low energy behaviour 
of the correlation functions. In the following we will call 
them the Luttinger parameters. They completely char- 
acterize the low energy physics. 

For excitations which do not change neither the num- 
ber of particles and the currents, the bosonic fields can 
be written in terms of boson creation and annihilation 
operators, bj, ^ and bv^k-, y = c,s 




+ (8) 



Jlvix) = i 



2LK„ 



k 



ikx 



-ikx 



(9) 



It is moreover possible to express the original fermionic 
fields in terms of the new boson fields by the so-called 
bosonization identity (see for instance [20,31,34,22]) 



a^o (27ra)V2 



(10) 



The fields (pr.a are recovered from Eq. (6). The (hermi- 
tian) operators F^^a are called Klein factors and obey the 
Clifford algebra {Fr.o-, ^r',^'} = 2(5r,r'<5(T,<T'- 

The same effective Hamiltonian, Eq. (7), can be ob- 
tained starting from the Hubbard Hamiltonian Eq. (3). 
In the optical lattice case, besides the back-scattering 
term, also an umklapp contribution could in principle be 
present, but it does not contribute away from half-filling. 



1. Luttinger parameters 

For exactly solvable models like the ID lattice model, 
the Luttinger parameters can be directly expressed in 
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terms of the microscopic parameters of the theory [35] . 
In a spin rotationaUy invariant Fermi-gas the quantity 
Ks = 1, so that the only independent parameters are 
Kp, and Vs,p [19,20]. 

One of the more peculiar consequences of the Hamil- 
tonian Eq. (7) is a complete separation of the dynamics 
of the spin and charge degrees of freedom. In general one 
has Uc ^ Us and thus the density (or charge) and spin 
oscillations propagate with different velocities. Only for 
a non-interacting gas, or for some accident, the two ve- 
locities are the same and equal to the Fermi velocity, 
Ui, = vp. This phenomenon is known as spin-charge sep- 
aration. Our aim is eventually to propose a way to see 
in a clear way such effect. 

We should mention that if gi± renormalizes to a finite 
value (in particular in a spin rotationally invariant sys- 
tem with bare gi± < 0, the parameter gets relevant by 
means of RG), one can have a gap in the spin sector, the 
charge excitations remaining massless [20,31]. 

In what follows we shall give some analytical results 
for the Luttinger parameters of systems described by the 
Hamiltonian Eq. (1) for a box and a harmonically trap- 
ping potential, with and without the underlying optical 
lattice. 

For weak interaction, we define the small parameter 
^ = g/nhvp. We get, up to the first order in ^, Uc = 
vp{l -\- ^/2), Us = vp{l - ^/2), Kc = l- C/2. The same 
results are valid in the presence of optical lattice, the only 
difference being in the definition of ^ = Ud/irhvp, with 
Fermi velocity vp = 2Jdsm{TTnd/2)/h [32]. 

When Vl.ct — 0, i.e. in the homogeneous case, the 
models described by the Eq. (1) and Eq. (3), are ex- 
actly solvable [36,37]. The solutions confirm that the 
low energy spectrum is Luttinger-like and so, at least in 
principle, it is possible to extract the Luttinger param- 
eters for any interaction strength. The exact solution 
[36,37] is suitable to study the strong interacting limit, 
i.e., ^ 3> 1. The parameters can be found using per- 
turbation theory around ^ ^ oo. Without the lattice 
one has: u^ = 2fnrn/3ml, = 1/2(1 -h 8 ln(2)/7r2^) 
and Uc = Trhn/m{l — 81n(2)/7r^5), where n = N/L is 
the density of the gas. When the lattice is present one 
obtains 



Jd 



Us = 



1- 



K. 



sin(27rn(i) 
7isin(7rnd/2)^ \ 2'Knd 

41n(2) sin(7rnd) \ 
TT^^ s\Ti{'imd/2) J 



Uc = 



2Jdsin{Tmd) 



1 - 



4 \n.{2)nd cos{TTnd) 
7rsin(7rn(i/2)^ 



(11) 

(12) 

(13) 



Some remarks can be made. First of all we should 

mention that the results without the lattice can be recov- 
ered from the previous expression expanding in the first 
(non-zero) order in low density n and substituting J — > 
h'^/2md'^, the new density being n = N/N^d N/L, 
with Ny, the number of sites. Indeed from the exact so- 
lution one gets that, in the limit of low density, the lattice 



case resembles the continuous case with a rcnormalized 
mass m* = {2Jd^/tT?)~^. Secondly, note that the results 
(put ^ ^ oo) confirm the intuition that when the repul- 
sion between the atoms of the two diS^erent species is very 
strong some properties of the gas are similar to those of 
an ideal single component gas of indistinguishable par- 
ticles. In particular the Fermi velocity is vp = irfin/m 
and the spin velocity goes to zero. Intuitively one can 
say that, somehow, the infinite (hard-core or on-site) re- 
pulsion plays the role of an effective Pauli principle for 
atoms in different internal state (i.e., different spin). The 
same is true for instance for strongly interacting bosons 
in ID, where there exist an exact mapping between hard 
core bosons and free fermions [38]. Note, however, that 
in our case the asymptotic value Kc = 1/2 [35], while 
is Kc = 1 for a non-interacting gas. It is well known 
that at zero temperature the Green's function of a Lut- 
tinger liquids present a peculiar power-law decay with 
non-universal, JsT-dependent exponents. At finite tem- 
perature the behaviour is still iiT-dependent [20]. This 
means that physical properties like compressibility, den- 
sity of states, momentum distribution clearly reveal the 
interacting nature of the gas. 



2. Inhomogeneous case: Local Density Approximation 

The inhomogeneous case will be studied using a local 
density approximation (LDA). In order to apply LDA 
to our trapped gas we should assume that the size R 
of the gas sample is much larger than the interparticle 
separation, i.e. R ^ kp^. Such condition is consistent 
with having a large number of particle, N ^ 1. The 
ground state of the system can be characterized using 
the Thomas-Fermi equilibrium condition; 



dEoin) 
dn 



IJ'-VL,a{x). 



(14) 



Where Eo{n) is the ground state internal energy of the 
gas per unit length, n the total density and /i is the chem- 
ical potential. The longitudinal potential in this case is 
Vl,(t{x) = l/2mw^a;^. The previous equation is just the 
expression of the fact that the energy cost of adding a 
particle to the system equals the chemical potential cor- 
rected by the local value of the external potential. In 
the Luttinger picture we assume that the spatial vari- 
ation of Ki, and Ui, originate only from the spatial de- 
pendence of the gas density, i.e., K^[n] — > K^[n{x)] and 
Uiy[n] Uv[n{x)]. 

The equation of motion for the fields conjugate fields 
representing the low energy excitations of the inhomo- 
geneous gas are obviously given by = [H, (j)^] and 
itli, = [H, H^] , where H is the Hamiltonian Eq. (7) with 
spatial dependent Luttinger paramters. Using the canon- 
ical commutation relation one can easily find: 



: K^{x)u^{x)Il„, 



d Uy{x) d 
dx K^{x) dx 



(15) 
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Let us first of all consider a non-interacting fermion 
gas, i.e. g = 0, without the optical lattice. This case is 
rather well-known (see for instance [14,15], where differ- 
ent approaches were used) . The density of the gas in the 
LDA in terms of the local value of Fermi momentum is 
n{x) = 2kF{x)/'!rh. The internal energy of the gas is just 
the density of the kinetic energy (a.k.a. quantum pres- 
sure) £^o(") = fi.^7r^n'^(a;)/24m. Substituting into Eq. (14) 
we find the ground stete Thomas-Fermi density profile 



riTFix) = no 4 



Rrpp 



(16) 



for \x\ < Rtf, and otherwise. We have indicated with 
no = {Sjim/h'^Tr^y^^ the central density of the trap and 
with Rtf = (2/x/ma;^)^/^ the Thomas-Fermi size of the 
cloud. The chemical potential /x = fkOxN/2 is fixed by 
the normalization requirement. 

For an ideal gas the equations of motion Eq. (15) can 
be solved. Indeed in this case the Luttinger parameters 
have the values u^, = vf, Kc = 1, where the Fermi ve- 
locity vf can be written in terms of the total density: 
Vf — h'Kn/2m. Substituting in the equation of motion 
Eq. (15) and using the density profile Eq. (16), we can 
write down the equations for the density and spin-density 
fluctuations (oc dx^i^)- 

-uIp. = ivl^il - - i'f'p., (17) 

where we defined the adimcnsional coordinate x = 
x/Rtf- The same equation has been derived in [15] us- 
ing the mean field equations of Kolomeisky et al. [39]. 
The solution is given by 



X^ Pvn 



= Al, sva.{(jJi,n/oJx arccosi) (18) 
-|- cos{ojvnl'^x arccosi). (19) 



By analyzing the boundary conditions one can find the 
discrete spectrum of eigenfrequencies. The result: w^n = 
UxTi [15] is the same both for the spin and the density 
modes. The first modes (n = 1) has eigenfunction p,^ ^ 
x/\/l — and it corresponds to harmonic oscillations 
of the total density and the total spin (dipole and spin- 
dipolc mode). 

For weak interaction the correction to the ground state 
energy can be obtained by averaging the interparticle in- 
teraction over the free ground state: 



£o = W-K^rfl2^m -|- gn^j^. 



(20) 



Then, in the spirit of LDA using the Eq. (14), one finds 

that the density, up to the first order in ^ = g/TThvF{0), 
uniformly decreases by dn{x) = —2gm/h^iT^. As ex- 
pected the (repulsive) interaction reduces the density. 
This simple conclusion breaks down close to the edges 
of the cloud x = ±Rtf , where the density goes to zero. 
Indeed we should require that g <C h'^n{x)/m, i.e., 



{Rtf - x)/Rtf < {gm/h^nof ~ 0(r)- (21) 

The Luttinger parameters can be calculated substituting 
the local Fermi velocity vf{x) = TTh{nTF{x) + Sn{x))/2m 
in the expressions obtained before for the homogeneous 
case with small interaction. In this way one finds 



_ wHriTFix) A _ Ac,sgm \ 
2m \ ■K'^ffiriTFix) J 



(22) 



where As = 3 and A,. = 1 and to the same order Kc = 
1 — gm/n'^h'^nTF{x). 

In the strong interacting limit from the exact solution 
one has 



Eo = r7r^n76m(l - Ah^nlTi{2)/mg). 



(23) 



The density profile for — > oo is still given by Eq. (16), 

but with no = n^o = (2m/i/7r^fi,^)^/^, /i = NfkUx and 
Rtf = Roo = (2/i/nza;^)^/^. The first order correction to 
the density is 5n{x) = 8h^ln{2)nTF{x)^/3mg. The spin 
velocity goes to zero as the interaction goes to infinity, 
the first order correction is simply 



Us{x) = TT h nTF{x)/3m g. 



(24) 



The density wave velocity and Kc have a non-zero value 
for infinite interaction and up to the first order correction 
one has Kc = 1/2(1 -|- 47i^ ln{2)nTF{x) / gm) and 



Uc{x) 



irhriTFix) 



m 



1 + 



4fi2 \n{2)nTFix) 
mg 



On the same lines we can obtain the results in presence 
of the optical lattice. For the non-interacting case {g = 
[/ = 0) the ground state energy of the gas is = 
2Jn — AJ suv{-Knd/2) /ttcI. From Eq. (14) one has that the 
density of the cloud is 



nTF = — ; arccos 
7ra 



1 



A(i 
2J ^ 



(25) 



for X = x/Rtf < 1 and elsewhere. From the previous 
expression we find the equations of motion for density 
and spin-density fluctuations in the ideal case 



(26) 



where wc defined fi{x) = ^/l -Ai/4J(l-52). For low 
density (i.e. /; 1) we have that the spectrum has 
the same structure as for the continuos case, b ut with 
a renormalized trapping frequency a;^,„ = nUxy^m/m*. 
This result is quite general and for instance the same cor- 
rection was found in [40] for a Bose-Einstein condensate 
in ID optical lattice. 

The first order correction, calculated averaging the in- 
teractionhoe[Coen term in Eq. (3) over the ground state is 
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5n = — /7nTF/27rJsin(7rn7'ird/2). The density and spin- 
density wave velocities up to the first order are given by 



Uc{s} = Vf{x) 



^ Ud / nTFdcos{TTnTFd/2) ^ 1 
2hvF \ sm{TrnTFd/2) tt 



(27) 



and to the same order the remaining parameter Kc has 
the same expression than in the homogeneous case but 
with a space dependent density given by the Eq. (25). 
In the strong interacting hniit the energy density of the 
ground state can be written as 



Eo 



2 J I n 

2 J2 ln(2) 
Ud 



sin{Trnd) 
Trd 



— sin(27rn(i) — (nd)^ 

TT 



(28) 



The density profile for J7 ^ oo is given by the Eq.(25), 
divided by 2. The rest of the calculation proceeds on the 
same line as above. We will use these results later to 
discuss spin-charge separation. 

We have shown how it is possible, under certain con- 
dition, to study as a Luttinger liquid an inhomogeneous 
ID Fermi ultra-cold gas. We gave some analytical results 
in the two limiting case of weak and strong interaction. 
This will be used in the next section to study a possi- 
ble way to see the equivalent of spin-charge separation in 
such system. 



IV. SPIN-CHARGE SEPARATION IN TRAPPED 
COLD-GASES 

A. Spectroscopy 

The separation of the spin and the charge modes in the 
spectrum suggests experimental ways to observe the spin 
charge separation in harmonically trapped atoms [27]. 
We will propose two related experiments to observe the 
separation of the spin and the charge degrees of freedom. 
One relies on the possibility to excite easily some normal 
modes of the trapped gas and to measure their frequency. 
The other one consist in really testing the dynamics of 
total density and relative density wave packets. Before 
discussing the experiments it should be noticed that the 
application of the Hamiltonian Eq. (7) to finite systems 
deserves special attention. As shown in [32] the bound- 
aries of a LL affect the physics of the lowest energy exci- 
tation e with eL/vF < 1. Setting L ~ Rtf we find that 
the excitations with LOun ~ ^iJx can have finite size correc- 
tions, which are not described by Hamiltonian Eq. (7) 
with the LDA values of the Luttinger parameters. Below 
we will consider LOj^n ^ (WKB approximation, see 
below) or small interactions. Hence we can neglect the 
finite size effect. Other finite size effect related to the TF 



approximation will appear (logarithmic correction), but 
they are not related to Luttinger physics. 

Assume we have the two different components of an 
interacting Fermi-system confined in external potentials, 
which act independently on each component. Then by 
modulating the two potentials simultaneously one in- 
duces oscillations of the total density. In turn, by slightly 
changing the confining potentials for the two components 
out of phase, only spin waves are excited. By observing 
the character and measuring the frequencies of the modes 
one can prove that the spin and the density modes are in 
fact the independent modes in the system and test the 
validity of the Luttinger description for such system. In 
FIG. 2 we show the numerical results of the level spac- 
ing for the spin and the charge modes for harmonically 
trapped atoms in the presence of optical lattice. Numer- 
ical results are based on the exact solution [37] together 
with the LDA. 
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FIG. 2. Level spacing (in units of Ux) of the spin (lower 
curves) and the charge (upper curves) modes for a harmoni- 
cally trapped Fermi gas in the presence of an optical lattice 
(see the text). The level spacing is shown as a function of 
^ = g/nhvF{0) and for different central densities: 
curves a: n(0) — 0.25/d (solid lines), 
curves b: n(0) = 0.42/ci (dot-dashed lines), 
curves c: n(0) = 0.58/d (dotted lines). 

For small interaction we can give an analytical expres- 
sion to the energy correction of the lowest modes for a 
trapped gas using the balance equation Eq. (14). The 
frequency of the (density) dipole mode is unchanged as it 
has to be, while the correction to the spin dipole mode di- 
verge logarithmically. The logarithmically divergent be- 
haviour is due to the localization of the excitations of 
a free gas close to the border (see before and [15]). To 
handle the divergence we use as cut-off the condition Eq. 
(21) and we find 



m^fiTFio) 



log 



-Kh^JlTFi^) 

2gm 
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The energy shift is negative and can be measured by com- 
paring the density and the spin osciUation of the atomic 
cloud. The logarithm in the previous expression tells 
us that in harmonic traps the requirement to apply the 
perturbation theory in ^ is stronger than in the homoge- 
neous case, namely it should be ^ log(l/^) <^ 1. We can 
say that the LDA is valid with logarithmic accuracy. 

The same behaviour for the first order correction to the 
spin-dipole frequency can be obtained in this regime also 
using sum rules. This was done for a 3D system in [16]. 
Generally speaking, it is possible to show, that an upper 
bound to the energy fuvp oi the lowest state excited by an 
operator F is given by the ratio h^cop = m^/mi, where 
rrife is the /cth moment of the dynamic structure factor. 
One can write the two moments as 



mi 



TO3 



\m[F\HUH,[H,F]m, 



(29) 



where H is the Hamiltonian and |0) the ground state 
of the system. We are interested in the excitations due 
to the spin dipole operator S = Xi — ^^j^ Xi and 
the Hamiltonian being Eq. (1). In the limit of small 
interaction one has 



j dx\dx 



n{x)\\ 



(30) 



Substituting for the ground state density n, the density 
profile Eq. (16) and using the cut-off Eq. (21), one finds, 
to logarithmic precision 



5ijJs = -to. 



16 



gm 



log 



gm 



(31) 



The Eq. (30) has nothing to do with Luttinger model 
and, hence, provides an alternatives derivation for the 
spin excitation frequency. The fact that the results agree 
with each other confirms both the applicability of the 
Luttinger model and the LDA. 

We will introduce now a WKB calculation of the en- 
ergy excitation, which we will allow us to calculate the 
energy shift of the spin excitations in the strong interact- 
ing regime. Suppose that the energy levels of a quantum 
system are charcterized by a certain quantum number n. 
For sufficently high quantum number the eigenfrequen- 
cies can be obtained from the WKB (Bohr-Sommerfeld) 
quantization condition. 



-xn 



p{x)dx = h'K{n + a), 



(32) 



where p{x) is the WKB momentum corresponding to a 
given energy, xq is the classical turning point and the 
constant a is the so called Maslov index. The accuracy 
of the WKB spectrum estimation is ~ l/(7rn)^. In our 
case we can fix a by comparing the WKB results and 



the exact solutions of Eqs. (15) for a weakly interacting 
gas. The spectrum of the excitations wc are considering 



is linear, e 



{x)p{x) , with the velocities given by the 



Eq. (22). Using the Eq. (32) we obtain the same sort of 
logarithmically diverging integrals as those we obtained 
in the perturbative calculation above. Once again to reg- 
ularize such diverging behavior we use the condition (21) 
and we find 



Cp^s = ^(« + 1) 1 - 



2gmAp 



log 



fi^riTF (0)7r 
mg 



(33) 



We find once more that the interaction dependent split- 
ting of the eigenencrgies is different for charge and spin. 

In the strong interaction limit, after integrating Eq. 
(32), we find that once again the density dipole mode is 
unchanged. The leading term of the energy of the spin 
dipole mode, which is zero in the limit ^ ^ oo, is given 

by 



S{V) 



hVn^{x)p{x) 
Sm'^g 



Substituting p{x) in the integral equality (32) and dealing 
with the logarithmic divergence as usual, we find, that 



huj{n + a) 



3gmlog{gm/h'^noo)' 



(34) 



where a ^ I. In the limit of the strong interaction the 
spin level spacing decreases a lot and is much smaller 
than that between the density waves (w^), as is clearly 
shown in FIG. 2. 



B. Wavepacket dynamics 

The difference between the spin and the charge veloci- 
ties allows to have spin and charge wave packets moving 
at diff'erent velocities. This can be seen in a sufficiently 
long trap. Localized wave packets of the spin and density 
excitations can be generated by sufficiently localized op- 
tical potentials of dimension £ such that R':$> £ ^ kp^, 
where R is the size of the atom cloud and the Fermi mo- 
mentum kp is essentially the interparticle distance FIG. 
3. 
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FIG. 3. Two component quasi-lD Fermi gas in a harmonic 
trap. A short laser pulse focused near the center of the trap 
can excites density end spin density wave packet. 

The procedure we have in mind is analogous to the 
MIT setup originally used to study propagation of sound 
waves in elongated condensates [41]. The motion of these 
wave packets is readily detected by internal state depen- 
dent laser probes. In FIG. 4 we show a simulate wave 
packets motion for the states \F = 1/2, Mp = ±1/2) 
of ^Li with interaction parameter ^ = 1. This figure 
was obtained solving numerically Eq. (15). The scat- 
tering length of these two states is magnetically tun- 
able and recently the tunability has been experimentally 
explored up to high magnetic field (see [42] references 
therein). The scattering length is positive for applied 
magnetic field B > 550G and presents a Feshbach reso- 
nance around 800G. We can, for instance, consider small 
scattering length of the order of a = 2QA. In order to 
have 1^ ~ 1, one can use N ~ 500 particles at trap fre- 
quencies ujx — IHz and uj± = 250kHz. On the other hand 
near the Feshbach resonance the scattering length can be 
larger by orders magnitude, with respect to the previous 
value. For a just ten times bigger than the previous one 
it is possible to have N = 1000 atoms at a trap frequency 
LUj_ — lOOkHz, leaving unchanged 
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FIG. 4. Wave packet dynamics for different times as func- 
tion of position (in units of the Thomas Fermi radius Rtf)'- 
spin-charge separation manifests itself in a spatial separation 
of the spin (solid line) and density (dashed line) wave packets 
(shown at half a trap oscillation period cot = 7r/2), which can 
be probed by a second short laser pulse at a later time. The 
parameters correspond to A'' = 10'^ ®Li atoms in a trap with 
uj = IHz (with coupling parameter ^ = 1, see text). 

An optical potential created by a far detuned laser 
light and acting in the same way on both the components 
amounts to an external potential ~ dx4'c and hence acts 
only on the density waves. An optical potential gener- 
ated by a laser tuned, e.g., between fine structure levels 
of excited Alkali states acts on the ground state spin in a 
way equivalent to an external magnetic field interacting 
with the spin density ~ dx4>s and thus act solely on the 
spin waves. One can also act at the same time coherently 
on both the charge and the spin waves. We analyze in 
some detail this case (the same can be done for the pre- 
vious cases). The idea is sketched in FIG. (5). We need 
to couple only one out of the two ("long-living") hyper- 
fine levels identified by |(t) = | t) J i); with another (in 
principle "short- living" ) state of the atom, |e). Let be 
the Rabi frequency which couples for instance | "f) and 
|e), S be the detuning, i.e. S = ujl — i-^^e, where uj^ is the 
laser frequency and fkuije is the energy separation between 
I t) and |e), and 7 be the spontaneous emission rate of 
the state |e). In the regime of large detuning, S ^ n,j 
the effect of the coupling is to introduce an energy shift 
= \n\^/S and a decay rate F = \^/S\^j. Thus, for 
not too long times {Tt <C 1), we have the Hamiltonian 
density Hoxt(a;) = {x)tp^^ {x)tp^ {x) acting on the gas. 
Eventually, if ri|(x) (i.e. the laser profile) is slowly vary- 
ing in space, to avoid scattering processes, q ~ 2kp, in 
terms of charge and spin fields one can write 

f dx 

Hcxt = J —^^T{x){dx(f>c{x) + dx(f>s{x)). (35) 

Clearly the previous Hamiltonian acts as localized (with 
the above remark) independent perturbations on both 
charge and spin densities. 
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FIG. 5. Internal level scheme of the atoms of the gas. The 
state I t) — \Mp = 1/2) is coupled via a laser to a (in general 
unstable) state |e) (left). The coupling being given by the 
Rabi frequency Q. On the right, we sketched, the effective 
action of the laser for large detuning 5 ^ $1,7 (see the text). 



V. RELAXATION PHENOMENA 

As stated before, the Hamiltonian Eq. (7) describes 
properly the fluctuations of a ID system only in the long 
wave q/kp <C 1, low energy regime. Higher order terms 
originating from, for example, the non-linearity of the 
fermionic spectrum appear which can mix the excita- 
tions. The first corrections are third power in Hi, and 
which lead to scattering of the excitations. This is 
precisely the case in the simple model described by Eq. 
(4). Introducing once again the operators Cr.k.cr, r — R,L 
and =t, i for the left- and right-movers in the momen- 
tum space, but this time keeping the terms up to the 
second order in k — kp , one has for the kinetic part 



r,k,a 



idle)k. 



{ark 



kpfclk^Ca.k,. 



(36) 



(37) 



We calculated the damping of the oscillations by the 
Fermi's Golden rule and the second quantized represen- 
tation for the phonon fields Eqs. (A2) and (A3). We 
considered the decay rate of an excitation with momen- 
tum q and energy ujq. In the limiting case ksT <^ hujq, 
a straightforward calculation gives 



T=0 



h^q^L/m^u. 



(39) 



This is the decay rate in a process where a particle 
with the energy hujq decays into a pair of particles with 
ll!i,lu2 < ujq. Note that ujq 3> LiJx, otherwise there are 
no final states for such decay instability (i.e. the lowest 
excitations are very stable at very low temperatures). In 
the case of fc^T > hcuq (but still T <^Tp) we find 



2„,3 



Tt q^Lk%T^/m^u 



(40) 



which corresponds to Landau damping, i.e., the contribu- 
tion of a process in which the damping occurs by scatter- 
ing a high-frequency excitation with hw'q ~ fc^T. Note 
that both results contain the size of the sample L. For 
sufficiently small temperatures (or high number of parti- 
cles) the excitations are only weakly damped. Indeed, for 
the lowest excitations Lmu/h ^ N ^ uoq ^ oj ^ ep/hN, 
thus substituting the previous expression for the decay 
rate we have 



T/iOq ^ [kBT/epf < 1. 



(41) 



In conclusion for sufhciently low tempertaure or large 
enough number of particles the excitations are only 
weakly damped by the third order correction to the other- 
wise quadratic Luttinger Hamiltonian. This means that 
in a degenerate quasi-lD Fermi gas the excitations can 
be observed for several trap periods. 



where a^i^i) — +(— ). The quadratic term in real space 
reads {dle)kpY.rM S ■ dxi>l,ai^)9xA,<Tix) and using 
the bosonization identity (see Appendix Eq. (Al)) one 
obtains 40r(9^e)A:j. J^r.a I dxid^tpa^a)^ ■ Introducing the 
charge and spin boson fields Eq. (6), we find that we have 
to add the following third order term to the Luttinger 
model density Hamiltonian 



Hint = ^M^mcj^cf 



V2Tr 

dx<t)jil + iicdx(t>cHc + n^a 

39.0c[(5,</.,)2 + n2]} 



'dIic{dx4>sHs + I\^sdx(t)s) 



2; 



(38) 



To study the relaxation phenomena we will consider, 
for simplicity, a single component gas (polarized system) . 
In the previous Hamiltonian we are left only with the 
term Hint = f?{dx<\>f' I ^mn. 

According to Haldane [19], the generic term due to 
the non-linearity of the spectrum is given by VJnt = 
7?i^/m / {dx4'Y, where 7 ~ 1. 



VI. SPINLESS LL WITH BOSONS 

In this section we will show how to realize a spinless LL 
Hamiltonian using "hard core" bosons in an optical lat- 
tice. The on site interaction of the bosons can be written 
as U{ni — l)n.i. Hence, in the limit of g — *■ cxd ([/ ^ 00) 
only the states rii = 0, 1 are allowed. We parameterize 
them as a = ± states of a particle with the spin 1/2. 
In terms of spin operators (Pauli matrices) the lattice 
(Hubbard) Hamiltonian can be put in the form 

H = jY.^4^i+i + h.c.) + vY,{a^r + i){at+^ + 1). 

i 

(42) 

Where V is due to the off-site interaction, which is 
present in the starting bosonic Hamiltonian. Usually one 
has V . In the easiest case, if d is the lattice period 
and 'w{x) the Wannier function, one has 



U = g dx\w(x)\'^ 



(43) 
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V = g j dx\w{x)w{x + d)\'^ . (44) 

The Hamiltonian Eq. (42) can be brought to a 
fermionic form with the help of Jordan- Wigner trans- 
formation: 

j<i 
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a- = 2a\a., - 1, 



(45) 



where the a^'s are fermionic operator. Substituting these 
expressions into the spin-chain Hamiltonian we get: 

B. = J^(a|ai+i +/I.C.) +4y^a^aiaJ_^iai+i. (46) 



Thus we have obtained a ID system of interacting spin 
polarized fermions starting with a system of bosons. 
Bosonizing the JW fermions we arrive at a LL Hamil- 
tonian similar to Eq. (7). 

We should also mention that, recently, was shown that 
it is possible to study typical fermionic correlation phe- 
nomena using 2-level bosonic atoms trapped in optical 
lattices [43]. 



VII. CONCLUSION 

In conclusion, we have described ID (Fermi) gases un- 
der different trapping configurations as Luttinger liquids 
by using the availbic exact solutions. We have shown 
how it is possible to gain some physical insights with this 
approach also in the non-interactig case (strictly speak- 
ing we cannot talk in this case of Luttinger liquids). In 
addition, we have proposed some experiments to test the 
validity of this approach and to see in a clear way the 
so-called "spin-charge separation" . 

We emphasize that a series of approximations has been 
made. We considered our system as an infinite one, or 
better finite with periodic boundary conditions, to ex- 
tract the Luttinger parameters. We commented on that 
for what concern an actual experiment. The finite-size 
correction can be, in any case with some effort, taken 
into account. The analysis of the inhomogeneous case 
was carried out within a local density approximation. For 
what wc gave the range of validity. 

In summary, cold quantum degenerate Fermi gases arc 
a promising avenue to observe Luttinger physics in a 
novel system in the laboratory. 
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APPENDIX A: BOSONIZATION IDENTITIES 

It is possible to represent the single fermion field oper- 
ators 'i/'r.o- {r = R,L and a =t, i) in terms of the boson 
field and U,y {v = c, s) (see for instance [20,34,22,31]) 
by the bosonization identity 



(27ra)i/2 



-ar»2\/7F<AB(I,),o 



(Al) 



where a/j(/,) = +(— )• The fields are recovered from 
Eq. (6). The (hermitian) operators Fr^a- are called Klein 
factors and obey the Clifford algebra {-Fr.o-, ^r',o-'} = 
26r,r'5a^a' . Wc remind that the Luttinger model is based 
on an infinite filled Dirac sea. The last correspondence 
must be handled with utmost care when a product of two 
(or more) fermion operators with the same argument are 
considcrc;d. In such cases one must implement the nor- 
mal ordering prescription, i.e., : A := A — (A)vac- For in- 
stance the right way to compute the densities is by using 
the limit prescription: p[x) = limg^o : '>p^{x)il){x + e) : 
[31]. 

It is useful to write the boson fields in terms of boson cre- 
ation and annihilation operators, bl j, and 




ikx I li „-ikx\ 



(A2) 



n.(x) = J ^ (^''.'^^''' - ^fc^"*'") 



L 



(A3) 



The eigenvalues of the operators iV^'s are the number 
of particles added to the ground state {y = c) and the 
"magnetization" (N^ —^\) of the system with respect to 
the ground state {y = s). The operators represent the 
respective currents. In the thermodynamic limit {L 
00 ) or for excitations which do not change neither the 
number of particles and the currents these operators can 
be neglected. Note that such "topological" excitations 
can be recognized also in the exact solutions of the models 
we discussed in this paper. 
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1. Backward and umklapp scatterng 

Using the previous identities it is possible to express 

also the backward scattering in terms of the bosonic fields 
and liu . In terms of fermionic operators the backward 
scattering is 



9l 2_y '-R,a\ 



(A4) 



Expressing the previous term in the real space and sub- 
stituting the fermionic fields with the Eq. (Al) one gets 

'Hi.ix) = igi±FR^Fii^FL^FLi cos(\/8.^,). (A5) 

The product of Klein factors can in this case substitute 
(with a kind of gauge choice) with +1 (see [20,31,34] for 
a comprehensive discussion). From the obtained expres- 
sion is quite clear that, when ^ij^ renormalizes to strong 
coupling, one has a gap for the spin excitations [44] . This 
is the case, e.g., for gi± < and Ks = 1. When the 
lattice is present also an umklapp term can be present. 
The "bosonized" form of such term is essentially like Eq. 
(A5), which (pc, instead of (ps in the argument of the co- 
sine function. But in this case a gap can be opened only 
at half filling [45]. 
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